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Abstract 

The apparent discovery of a Higgs boson with mass ~ 125 GeV has had a significant impact on the 
constrained minimal supersymmetric extension of the Standard Model in which the scalar masses, 
gaugino masses and tri-linear A-terms are assumed to be universal at the GUT scale (the CMSSM). 
Much of the low-mass parameter space in the CMSSM has been excluded by supersymmetric 
particle searches at the LHC as well as by the Higgs mass measurement and the emergent signal 
for B s — > Here, we consider the impact of these recent LHC results on several variants 

of the CMSSM with a primary focus on obtaining a Higgs mass of ~ 125 GeV. In particular, we 
consider the one- and two-parameter extensions of the CMSSM with one or both of the Higgs masses 
set independently of the common sfermion mass, mo (the NUHM1,2). We also consider the one- 
parameter extension of the CMSSM in which the input universality scale Mj n is below the GUT scale 
(the sub-GUT CMSSM). We find that when M{ n < Mqut large regions of parameter space open up 
where the relic density of neutralinos can successfully account for dark matter with a Higgs boson 
mass ~ 125 GeV. In some of these regions essential roles are played by coannihilation processes that 
are usually less important in the CMSSM with Mj„ = Mqut- Finally, we reconsider mSUGRA 
models with sub-GUT universality, which have the same number of parameters as the CMSSM. 
Here too, we find phenomenologically viable regions of parameter space, which are essentially non- 
existent in GUT-scale mSUGRA models. Interestingly, we find that the preferred range of the 
^4-term straddles that predicted by the simplest Polonyi model. 
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1 Introduction 



Supersymmetry is one of the best-motivated extensions of the Standard Model, including 
the naturalness of the electroweak scale, the appearance of a candidate for the cold dark 
matter if i?-parity is conserved, unification of the fundamental interactions, the prediction 
of a Higgs boson weighing ^ 130 GeV, and the stabilization of the electroweak vacuum if 
the Higgs boson weighs ^ 127 GeV, as seems now to be the case. For these and other, 
more theoretical reasons, much attention has been focused on models with supersymmetric 
particles weighing ^ 1 TeV, and notably the simplest, minimal supersymmetric extension 
of the Standard Model (the MSSM). In particular, many studies have been made of the 
simplified variant in which the soft supersymmetry-breaking scalar masses, gaugino masses 
and tri-linear A-terms are assumed to originate at some high mass scale and be universal at 
the supersymmetric GUT scale (the CMSSM) [lfj5]. 

However, there is as yet no direct evidence for supersymmetry. Some years ago, hopes 
were high that the deviation of the experimental measurement of the anomalous magnetic 
moment of the muon, — 2 [6], from the prediction of the Standard Model might herald 
the appearance of low-mass supersymmetric particles accessible to the initial run of the 
LHC However, this has not yet transpired, and the continuing absence of supersymmetric 
particles at the LHC [7| is putting increasing pressure on the CMSSM and rendering the 
supersymmetric interpretation of g^ — 2 somewhat dubious in this and related models. 

The ATLAS and CMS experiments have recently discovered a new boson with spin ^ 1 
and couplings similar to those of the Standard Model Higgs boson |8]. This is compatible 
with such simple supersymmetric models, which predict not only the existence of a Higgs 
boson with mass ^ 130 GeV, but also that its couplings to Standard Model particles should 
be indistinguishable from those in the Standard Model at the present level of accuracy. 
Using FeynHiggs [9] there is an uncertainty ~ ±1.5 GeV in the prediction for the Higgs 
mass for any given set of input CMSSM parameters, but the measurement of a mass ~ 125 
to 126 GeV does disfavour relatively low sparticle masses. Taken together with the absence 
so far of supersymmetric particles at the LHC, the Higgs mass measurement excludes a large 



range of low-mass CMSSM parameters 10 



14 



In parallel, increasing pressure is being placed on models with large tan (3 by the strength- 



ening experimental constraint on the decay B s — >■ [i ,ft~ 15 , in particular the recent de- 
tection by LHCb of B s — > with a value close to the Standard Model prediction [16]. 
Currently, the 95% upper limit on the branching fraction relative to the Standard Model 
value is 1.65 from LHCb alone, improving to 1.50 when combined with other experiments. 

These developments motivate the study of generalizations of the CMSSM, the search 
for viable alternatives, and the exploration of their observational signatures, with a view to 
future searches at the LHC and elsewhere. 

In this spirit, here we consider various one- and two-parameter extensions of the CMSSM. 
One option that has often been studied is to allow one or both of the soft supersymmetry- 
breaking contributions to Higgs scalar masses to be independent of the common sfermion 
mass, mo (the NUHM1,2) [P^O]. These offer some extra possibilities for accommodating 
the mass of the Higgs boson and the absence (so far) of supersymmetric particles at the LHC 
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while yielding the appropriate relic cold dark matter density. 

Another way to extend the CMSSM with one additional parameter is to allow the input 
universality scale to be lower than the GUT scale [21] (the sub-GUT CMSSM). In this case, 
the sparticle spectrum is compressed relative to that in the CMSSM, because of the reduced 
amount of renormalization of the soft supersymmetry-breaking masses if they are universal 
below the GUT scale. This in turn implies that relatively large ranges of sparticle masses 
are compatible simultaneously with a Higgs boson mass ~ 125 GeV and an appropriate relic 
dark matter density, often involving novel coannihilation mechanisms not relevant in the 
CMSSM. 

As a corollary, one may impose an additional constraint on the sub-GUT extension of the 
CMSSM and retain many of its phenomenological advantages. In particular, one may impose 
the specific relation Aq = Bq + mo between the trilinear and bilinear soft supersymmetry- 
breaking parameters A and B that is characteristic of minimal supergravity (mSUGRA) 



theories 22 -24 . Intriguingly, we find acceptable models for a range of A values that includes 
the specific value Aq = (3 — \/3) x 7713/2 found in the original Polonyi model [25]. 

A generic feature of the NUHM1,2 and such sub-GUT models is the increased importance 
of coannihilation effects, that are only important along narrow strips in the CMSSM param- 
eter space. Accordingly, in this paper we have paid particular attention to the treatment of 
coannihilation processes, using a more complete treatment than in our previous analyses of 



theNUHMl,2 18,20 and sub-GUT models 21 



The layout of this paper is as follows. In Section 2 we recall relevant features of gravity- 
mediated models such as mSUGRA and the CMSSM, as well as pointing out some important 
features of their extensions to the NUHM and sub-GUT models. In Section 3, we discuss 
in more detail the impact of the LHC constraints and briefly review the current constraints 
on the CMSSM and mSUGRA. Section 4 discusses the impact of the LHC constraints on 
the NUHM1,2, Section 5 discusses sub-GUT generalizations of the CMSSM, and Section 6 
discusses the specific case of sub-GUT mSUGRA. Finally, Section 7 summarizes our principal 
conclusions. 



2 mSUGRA, the CMSSM and Generalizations 



2.1 mSUGRA 

We start by recalling relevant aspects of minimal supergravity (mSUGRA) models, which 



have a flat Kahler potential and a low-energy scalar potential that can be written as 22 - 24 

2 



V 



dW 



+ (A W {3) + B W {2) + h.c.) + mj /2 cj) 



i i 



where W is the superpotential for matter fields, 

W = (ye^Le + y d E x Qd c + y u H 2 Qu c ) + f iH 1 H 2 . 



(1) 



(2) 



and the SU(2) indices have been suppressed. In ([l|, is the trilinear part of the su- 

perpotential, is the bilinear part, and m 3 / 2 is the gravitino mass. This is the order 
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parameter for local supersymmetry breaking, and in mSUGRA one finds scalar mass uni- 
versality, m = 777.3/2- This condition is applicable at some input renormalization scale, M in , 
which is usually associated with the grand unification scale. In addition, in mSUGRA there is 
a relation between the tri- and bilinear supersymmetry breaking terms, A Q = B + m . With 
a minimal choice for the gauge kinetic function, gaugino mass universality characterized by 
7774/2 is also obtained, which is usually imposed at the same input scale M in . 

Minimization of the Higgs potential leads to two vacuum conditions at the weak scale, 
which can be expressed as 



m\ - m\ tan 2 /3 + kmUl - tan 2 0) + A^ 

V = o - : — vj2) > ( 3 ) 



.2 _ 

tan 2 /? - 1 + A^ 



and 

B(i = --{m\ + m\ + 2/i 2 ) sin 2/3 + A B , (4) 

where mi i2 are the soft supersymmetry-breaking Higgs masses (evaluated at the weak scale), 
tan/3 is the ratio of the two Higgs vacuum expectation values, and Ab and A^' are loop 
corrections 26 . As a result, an mSUGRA model can be defined in terms of just 3 continuous 
parameters: mi/ 2 , m and A , and the sign of \x. Together, these determine tan/3 and the 
magnitude of \x via the vacuum conditions (|3j El). 

2.2 The CMSSM 

The CMSSM is the most thoroughly studied of all classes of constrained supersymmetric 
models [lJ-jH]. It is a two-parameter generalization of mSUGRA, in which the relation between 
A and Bq is dropped, allowing one to treat tan/3 as a free parameter, as well as dropping 
the relationship between the universal sfermion mass mo and the gravitino mass, which is 
possible if the supergravity Kahler potential includes non-minimal kinetic terms. Thus, the 
CMSSM is defined by the sign of /i, 777,1/2, m o, A , tan/3, and 7773/2, though the latter is 
usually ignored, e.g., because it is assumed implicitly to be irrelevantly large. We note 



that the CMSSM can be directly related to mSUGRA 27 through an extension in which 



terms proportional to are added to the Kahler potential as in the Giudice-Masiero 



mechanism 28 . For example, a contribution to K of the form 

AK = c„H x Hz + h.c. , (5) 

where ch is a constant, and i?i j2 are the usual MSSM Higgs doublets, affects the boundary 
conditions for both \i and the B term at the supersymmetry breaking input scale, M in , in a 
manner similar to the CMSSM. 



2.3 The NUHM1,2 



One of the possible one-parameter extensions of the CMSSM is the oft-studied NUHM1 17 



18 , in which the soft supersymmetry-breaking Higgs masses are taken to be equal to each 
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other, but are allowed to differ from the otherwise universal sfermion mass, ni\ = ni2 7^ mo- 
We also discuss here a two-parameter extension of the CMSSM, namely the NUHM2 [20] in 
which TOi 7^ m 2 7^ m , in general. 



2.4 Sub-GUT versions of CMSSM and mSUGRA 

Another one-parameter generalization of the CMSSM is provided by models in which the 
input scale for supersymmetry breaking universality differs from the GUT scale. Although 
M in might be chosen above the GUT scale |27j|29}|30] , such a choice would introduce many 
more parameters associated with different GUT models, and its study lies beyond the scope 



of this work. Here we concentrate on 'sub-GUT' models in which Mi n < Mqut 21 , which 



can be related to mirage unification models 31 



Both the NUHM1 and sub-GUT CMSSM are six-parameter models (five if one neglects 
the gravitino mass, e.g., by assuming that it is irrelevantly large). As we shall see, both 
the NUHM1,2 and sub-GUT models allow a suitable relic dark matter density and a Higgs 
mass of 125 GeV to be achieved simultaneously. This is also the case in sub-GUT versions of 
mSUGRA. This is a four-parameter class of models (m^, ^3/2, A , and M in ). Interestingly, 



we find that such models are phenomenologically viable in a re 
A , which straddles the Polonyi value A = (3 — y/3) x m 3 / 2 



24 



atively restricted range for 
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As an appetizer for the subsequent discussion, we display in the left panel of Fig. [T] 
the evolution of the sparticle spectrum with Mj n < Mqvt for a typical Polonyi mSUGRA 
scenario with m\M = mo = 2000 GeV. At these values of (mu2, mo), tan/3 varies from ~ 18 
when M in = Mqut to ~ 43 when M in = 3 x 10 s GeV: we do not find consistent electroweak 
vacuum solutions much below this value of M in , as indicated by the mauve shading in Fig. [T] 
As M in decreases, we see that the spectrum compresses until M in ~ 10 10 ' 5 GeV with, in 
particular, the approach of m X3 and m X2 to m x . This has the general effect of enhancing 
the importance of coannihilation processes that are less suppressed by Boltzmann factors 
~ exp(— Am/T). We also note that m x ~ m^/2 when M in ~ 10 12 GeV, indicating that 
the parameters traverse a rapid s-channel annihilation funnel region. There is level crossing 
in the neutralino spectrum when M in ~ 10 10 ' 5 GeV, with the composition of the lightest 
supersymmetric particle (LSP) x making a transition from almost pure bino to almost pure 
Higgsino, as seen in the right panel of Fig. [TJ For M; n ^ 10 10 5 GeV, the LSP and the next- 
to-lightest supersymmetric particle (NLSP) are nearly degenerate Higgsinos, with important 
implications for annihilation rates. We also note that their masses approach m^/2 again 
just below 10 9 GeV, so that rapid (co) annihilations via s-channel H/A poles again become 
important. 

This example illustrates just some of the variety of (co) annihilation phenomena that can 
control the relic LSP density in sub-GUT models. More of the possible richness will be seen 
in subsequent figures. 



4 




11 12 13 14 15 16 ' 8 9 1U 11 " 1J 

log 10 (M in /GeV) log 10 (M in /GeV) 



Figure 1: The evolution of (left) the sparticle spectrum and (right) the composition of the 
LSP x as functions of M in in a specific sub-GUT Polonyi mSUGRA scenario with muz = 
m = 2000 GeV. 

3 The Impacts of Recent LHC Results on the CMSSM 
and mSUGRA 



The following are the principal experimental, phenomenological and cosmological constraints 
on supersymmetric models that we consider in this paper: 



The preferred cosmological cold dark matter density range is 32 : QcDhih = 0.112 ± 
0.006, where h is the present Hubble expansion rate in units of 100 km/s/Mpc. We 
emphasize models in which the lightest supersymmetric particle (LSP), assumed here 



to be the lightest neutralino, \i is the dominant source of cold dark matter 33 . In 
general, supersymmetric models in which the LSP density is not excessive are squeezed 
closer to the boundaries where the LSP becomes charged or there is no consistent 
solution to the electroweak vacuum conditions ([3j [4]), and do not have very different 
values of the input parameters. 

The strongest upper limit on the spin-independent cold dark matter scattering cross 



section is from XENON100 34 , which excludes the focus-point region of the CMSSM 



and some analogous regions in models with more parameters |10] . 

The anomalous magnetic moment of the muon, — 2, will be treated as an optional 
constraint j6], unlike the previous constraints, which are considered to be mandatory. 
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Constraints are provided by flavour physics, in particular by b — > 57 [35] and B s 



— >■ 



15 16 , the latter now being dominated by the measurement by LHCb and upper 



limits from the other LHC experiments. 

ATLAS and CMS provide limits from searches for missing-energy events at the LHC, 
in particular the ATLAS search with ~ 5/fb of data at 8 TeV [7 36 , which is currently 
the most sensitive public result. 



The measured mass of the (lightest) Higgs boson is taken as Mh = 125.7 ± 1.0 GeV 
We recall that the theoretical calculation uncertainty of Mh for any given set of input 
parameters for the CMSSM or a similar model is typically ~ 1.5 GeV, so any parameter 
set yielding Mh > 124 GeV could probably be regarded as acceptable. 



3.1 Impacts on CMSSM models 

Fig. [2] displays the interplays of these constraints in the (777,1/2, tuq) planes of the CMSSM for 
four choices of the input parameters Aq and tan /3, assuming in all cases that fi > The 
left panels of Fig. [2] are for tan/3 = 10, and the right panels are for tan/3 = 40. The upper 
panels are for Aq = 0, and the lower panels are for Aq = 2.5 mo- In each panel, the region 
at high m.1/2 an d low m where the 77 is the LSP is shaded brown, and the region at low 
777.1/2 an d high tuq where there is no consistent solution to the vacuum conditions ([3j |4j) is 
shaded mauve. The regions excluded by b — > 57 are shaded green, those favoured by — 2 
are shaded pink, and those favoured by Q x h 2 are shaded dark blue ^ The LEP chargino 
exclusion is shown as a near- vertical dashed black line at small m^j |37|. 

Turning to the LHC constraints, contours of Mh are shown in all panels as red dash- 
dotted lines, and the ATLAS exclusion from the search for events with missing transverse 
energy (MET) is shown as a solid purple line [^] We also display as solid green lines three 
contours of BR(B S — > /i + fi~)/BR(B s — > [i + [i~)sM = 1.65,1.5 (the present 95% CL upper 
limit from LHCb and combined experiments), 1.3 (the 68% upper limit from LHCb). For 
reference, where relevant we also show the contours for 1.0 (the Standard Model value) and 
0.85 (corresponding to the central value seen at LHCb). The same shadings and line styles 
are used in all subsequent figures. 

Concerning the impacts of the LHC constraints, we see in the upper panels that Mh is 
in general too small if Aq = 0. We also see that the MET constraint allows only the upper 



1 Historically, this sign of fj, has been favoured by the supersymmetric interpretation of — 2 (we display 
g^ — 2-compatible regions in Fig. [2] but not in subsequent figures) and to facilitate compatibility with b — > sj. 
However, since the LHC constraints now disfavour regions of the CMSSM parameter space that 'explain' the 
g^ — 2 discrepancy (as seen in Fig. [2j), this assumption should be taken with a grain of salt. 

2 Because of the combination of recent experimental constraints, we are forced to consider substantially 
larger ranges in the displayed supersymmetric parameters such as m and Therefore, here and in 

many other figures, for reasons of visibility we shade wider strips where 0.06 < Q x h 2 < 0.2. 



3 As discussed in 10 , this limit is generally stronger than other direct searches for supersymmetry at the 



LHC in the models studied, and is largely independent of tan/3 and Aq. Based on the analysis in 39 , in 
the upper left panel of Fig. [2] we have assumed that the constraint published in 36 is independent of mi/2 
when extrapolated to lower m . 
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tan p = 10,A = 0,(X>0 tan (3 = 40, A = 0,|x > 




m 1/2 (GeV) m 1/2 (GeV) 



Figure 2: The CMSSM 
A = (upper) and A Q = 



(wi/2) Wo) planes for fi > 0, with tan/3 = 10 (7e/^ and 40 (right), 
2.5 m (lower), as calculated for m t = 173.2 GeV using the latest 



version of the SSARD code '381 . The interpretations of the shadings and contour colours are 
described in the text. 



end of the co-annihilation strip close to the f\ LSP boundary Q and is incompatible with a 
supersymmetric resolution of the g M — 2 discrepancy. The B s — > fi + \i~ constraint has no 

4 We recall that the focus-point strip is excluded by the XENON100 upper limit on spin-independent dark 
matter scattering. 



7 



impact for tan (3 = 10, but for tan (3 = 40, Aq = its restriction on the coannihilation strip 
is stronger than that of the MET constraint. 

In the lower panels, we see that is generally larger when A = 2.5 m . However, for 
tan f3 = 10, Mh does not grow above ~ 121 GeV, whereas for tan (3 = 40, Aq = 2.5 wiq there is 
compatibility for rrii/2 ?t 1 TeV along the coannihilation strip. This region is also compatible 
with the LHC MET constraint, but not with the supersymmetric interpretation of g M — 2. 
Note that for the larger values of A , we see another brown shaded region in the upper left 
where the stop becomes the LSP (or tachyonic). To the right of this boundary there is a 
stop coannihilation strip which occurs at relatively low when tan (3 = 10 and is in the 
region excluded by b — > 37 when tan (3 = 40. Once again, the B s — > fi + fi~ constraint has no 
impact for tan/3 = 10, and requires large mi/2 for tan/3 = 40. The disconnect between the 
LHC results and the potential solution to the — 2 discrepancy is severe enough that we 
will not show the — 2 compatible regions in future plots. 



A more complete survey of CMSSM models with other values of A is given in 14 , but 
this brief summary serves to illustrate how the LHC constraints prefer larger values of tan j3, 
Aq and m 1 / 2 . 

3.2 Impacts on mSUGRA models 



In Fig. |3jwe show a pair of mSUGRA (mi/ 2 , m ) planes for /i > 0. We recall that, there is just 
one free parameter in addition to these, which may be taken as Aq: the electroweak vacuum 



conditions (|3j |4J) may be used to determine tan (3 at each point on the plane 40 . In the left 
panel, we adopt the value Aq = (3 — m found in the original Polonyi model. Contours 
of fixed tan f3 are shown by the grey solid curves in increments of 5, as labelled. Also shown, 
as a diagonal and solid light blue line, is the contour where 777,3/2 = min(m x , m^), below 
which the gravitino is the LSP. Another diagonal line (brown dotted) shows the contour 
where the lightest neutralino mass m x is equal to the mass of the lighter stau, m^. The 
latter appears below the gravitino LSP line, so the fj is never the LSP in this particular 
model. As a consequence, only the dark blue shaded region at low m 1( / 2 above the light 
blue line corresponds to neutralino dark matter with the WMAP density. The dark blue 
shaded region below the light blue line corresponds to the gravitino LSP having the WMAP 
density it would have if there is no non-thermal contribution to the gravitino density. Here, 
the gravitino density is determined from the (co) annihilations of the NLSP - either the 
neutralino or (if below the dotted line) the lighter stau - followed by decay of the NLSP 
into the gravitino LSP, resulting in Q 3 / 2 h 2 = (m 3 / 2 /m x/?1 )f2 Xi f 1 /?, 2 . We recall that there are 
important cosmological and astrophysical constraints on long-lived NLSP decays, which we 
do not discuss here. 

Concerning the LHC constraints, we note that Mh is always far below the LHC measure- 
ment, which is a much stronger constraint than that from the MET search. Both constraints 
are incompatible with the supersymmetric interpretation of — 2. Since tan (3 < 35 in the 
portion of the mSUGRA (7771/2, m ) plane for Aq = (3 — \/3) uiq displayed in the left panel 
of Fig. [3j the B s — > constraint has no impact. 

The right panel of Fig. [3] shows a similar analysis for the case Aq = 2rriQ. Comparing 
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100 1000 2000 100 1000 2000 

m 1/2 (GeV) m 1/2 (GeV) 

Figure 3: The mSUGRA (mi/ 2 ,mo) planes for /x > and Aq = (3 — \/3)m (left), A = 2mo 
(right). The interpretations of the shading and contour colours are described in the text. 

with the left panel, we see that the values of tan j3 determined by the electroweak vacuum 
conditions are generally larger. For this reason, there is a region at low mi/2 that is excluded 
by b —>• S7, and the B s — >■ constraint rules out a significant portion of the (m^, mo) 

plane. We also see that the LHC Mh constraint is respected for mo ^ 1000 GeV. In this case, 
there is a region where the stau is lighter than the lightest neutralino and both are lighter 
than the gravitino: this f\ LSP region is shaded dark brown. In the lower right wedge at 
large mj/j and small mo, below the brown excluded region, the gravitino is once again the 
LSP and the f\ is the NLSP. In this case, above the f\ LSP region, we do find a WMAP 
co-annihilation strip, which extends to m^j ~ 1100 GeV. However, this strip has Mh ^ 123 
GeV and fails to respect the B s — > fi + n~ constraint. 

This summary serves to illustrate a broader incompatibility between the LHC constraints 
and mSUGRA, traceable to the extra constraints in this model compared to the CMSSM. 
We note that a more complete statistical analysis which included early LHC results already 



showed incompatibilities in the mSUGRA constructions 41 



4 The NUHM 

4.1 Results for NUHM1 Models 

We recall that the extra parameter introduced in the NUHM1 allows either itla or /i to be 
treated as a free parameter when solving the electroweak vacuum conditions (|3j |1J) for any 
set of fixed values of m 1 / 2 ,m , A and tan/3. Fig. Ill displays some representative (ffli/2, 
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planes for fixed Aq, tan /3 and tua, and Fig. [5] displays analogous planes for fixed values of /i. 



tan B = 10, m A = 1000 GeV, A = 2.5 m , u. > tan B = 30, m A = 1000 GeV, A n = 2.5 m n , u. > 

3000 n ■ i p- £ 1 1 ' — ^ 1 V-J- 1 3000-. A 



> 

O 



2000 



1000 



m h = 124 GeV 




m 1/2 (GeV) 



> 

o 



2000 



1000 



1500 




m 1/2 (GeV) 



1500 



Figure 4: The NUHM1 (rni/2,Wo) planes for fi > 0, itla = 1000 GeV, Aq = 2.5 mo and 
tan/3 = 10 (left) and = 30 (right). The interpretations of the shading and contour colours 
are described in the text. 



When rriA is held fixed as in Fig. [4], one finds an s-channel rapid annihilation funnel at a 
value of mi/2 slightly greater than itla (we recall that the funnel appears when 2m x ~ tua, 
and that m x ~ 0.43mi/ 2 )- This is seen as a vertical structure at mi/ 2 S> 1100 GeV H that is 
comfortably compatible with the LHC MET searches. For tan j3 = 10, the funnel does not 
extend past Mh ~ 123 GeV, whereas it extends almost up to Mh = 124 GeV for tan = 30. 
In both cases, the Higgs mass is probably too small to be compatible with the LHC result. 
Again, on the other hand, since the funnels shown here were chosen to occur at relatively 
large mw 2) there is no conflict with the LHC MET searches. For tan/3 = 30, the branching 
ratio of B s — > fi + n~ is somewhat higher than the Standard Model value, but it is still within 
the current 95% CL upper limit. At tan/3 = 40 (not displayed), there is increased tension 
with B s — )■ ji + n~ as the funnel sits very close to the 95% CL upper limit (a ratio of 1.5). In 
addition, at the higher value of tan (3, the end point of the funnel rises up only to 123 GeV. 

Large regions of acceptable parameter space are found in the (m 1 / 2 ,m ) plane when \i 
is held fixed, as seen in Fig. [5] where fi = 500 GeV. Here we have again fixed tan = 10 
and A = 2.5mo in the left panel, but have taken tan (3 = 40 and Aq = 2m in the right 
panel. For larger Aq/uiq, the fi LSP region (shaded brown) collides with the region (shaded 
mauve) where there is no radiative electroweak symmetry breaking. At low mi/2, there is a 
fx coannihilation strip which has low Mh when tan /3 = 10, and is for the most part in conflict 
with b — > sj when tan/3 = 40. At larger m^, there is an increasing Higgsino component in 

5 The width of the structure is enhanced for visibility by colouring the area where 0.06 < il x h 2 < 0.2. 
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tan [3 = 10, jj, = 500 GeV, A = 2.5 m tan p = 40, = 500 GeV, A = 2.0 




100 1000 2000 100 1000 2000 



m 1/2 (GeV) m 1/2 (GeV) 

Figure 5: The NUHM1 (mi/2, m o) planes for fi > 0, fi = 500 GeV, tan/3 = 10 and 

A = 2.5 m (left) and tan/3 = 40 and A Q = 2.0 m (right). The interpretations of the 
shading and contour colours are described in the text. 



the LSP and, as result, an increased annihilation cross section. This results in the relatively 
thick vertical strips seen at ~ 1200 GeV in both panels, again comfortably consistent 
with the LHC MET searches, where the transition towards a Higgsino-like LSP causes the 
relic density to fall into the range favoured by WMAP. At still larger mi/2, the relic x density 
is too small, but this region is not excluded if there is another dark matter candidate. 

For the fixed parameter choices in Fig. [5j this transition strip lies beyond the current 
LHC MET search. In this case, it is possible to obtain a Higgs mass in excess of 125 GeV 
even for tan /3 = 10, and much of the transition strip lies in excess of = 126 when 
tan/3 = 40. As seen from the contours showing the branching ratio of B s — > fi + fi~, this 
region is also compatible with the recent LHCb result. 

There are of course many possible slices that one can make through the NUHM1 param- 
eter space, and we only show a few of the more interesting ones here. In Fig. [6j we show 
examples of (/i, m^) planes with fixed m = 1000 GeV and A = 2.5 m , tan/3 = 10 on 
the left and tan/3 = 30 on the right. In both panels there is a brown shaded region at low 
mi/2 and small \fi\ corresponding to a stop LSP (or even a tachyonic stop). Just above the 
shaded region, there is a thin strip where the relic density falls in the WMAP range due 
to stop co-annihilations. However, if it is not excluded by b — > s'j, the value of is too 
low. Jutting out from this region, are two antenna-like strips with m\/2 ^ 800 GeV and 
hence compatible with the LHC MET constraint, which correspond again to the transition 
between bino and Higgsino-like LSPs we met in the previous figure. In the region between 
the two antennae, the relic density is too small. For tan (3 = 10, Mh is in excess of 124 GeV, 
and for tan /3 = 30, it is in excess of 125 GeV. For tan /3 = 30, the brown shaded regions at 
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large mi/ 2 depict regions with a f\ LSP. At still higher tan/3, these would descend and cover 
much of the plane. 




\i (GeV) |j, (GeV) 



Figure 6: The NUHM1 (//, mi/2) planes for mo = 1000 GeV, Aq = 2.5 mo, and tan/3 = 
10 (left) and tan/3 = 30 (right). The interpretations of the shading and contour colours are 
described in the text. 

It is interesting to note that in the right panel of Fig. [6] the branching ratio of B s — > 
falls below the Standard Model value, particularly when fi < 0. Here we find a region with 
acceptable relic LSP density, a Higgs mass M h ~ 125.5 GeV, acceptable b — > sj, and 
B s —> fi + fi~ at 0.9 its Standard Model value, close to the central value measured by LHCb. 

In Fig. [7] we show one more example of NUHM1 planes, a pair of (//, mo) planes with 
fixed m.1/2 = 1000 GeV (and hence LHC MET-compatible) and Aq = 2.5 mo, with tan f3 = 10 
on the left and tan/3 = 30 on the right. We again find regions with either a ti LSP (high 
mo) or a f\ LSP (low m ) with corresponding adjacent ti or f\ coannihilation strips. These 
coannihilation strips have either a low Higgs mass or are in conflict with b — > 57. At larger 
values of we also see regions where sleptons of the first two generations can become 
the LSP (coloured black). In both panels, there are vertical transition strips with suitably 
large values of Mh- The constraint from b — > 57 becomes particularly strong for fi < and 
tan/3 = 30, while the transition strip for /i > has a branching fraction for B s — > fi + fi~ 
close to the Standard Model prediction. 

These examples illustrate that the NUHM1 has the flexibility to accommodate the LHC 
constraints in a quite generic way, exploiting the extra flexibility in the Higgs sector of this 
model. 
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-2000 -1000 1000 2000 -2000 -1000 1000 2000 



\i (GeV) \i (GeV) 

Figure 7: The NUHM1 (/i, m ) planes for m 1( / 2 = 1000 GeV, A = 2.5m , and tan/3 = 
10 (left) and tan/3 = 30 (right). The interpretations of the shading and contour colours are 
described in the text. 

4.2 Results for NUHM2 Models 

There are still more planes that can be constructed in the NUHM2, with its extra free 
parameter. Here we show only one pair of (/i, ttia) planes that are LHC MET-compatible, 
with fixed mi/2 = rn = 1000 GeV and A = 2.5mo and, as in previous planes, tan/3 = 10 
on the left and tan /3 = 30 on the right. In both these particular examples, we see cruciform 
regions of acceptable relic density formed by the combination of a funnel, around fixed 
rriA — 900 GeV, and a transition to a Higgsino LSP when \fj,\ < 500 GeV. The constraint from 
b — > 57 disfavours more of the parameter space with /i < 0, particularly when tan = 30. For 
this value of tan /3, the constraint from B s — > excludes the surviving portions of the relic 

density strips below rj 1000 GeV when \x > 0. However, at large m^, BR(5 S — » 
is sufficiently small and ~ 125 GeV. In the case of tan /3 = 10, ^ 124 GeV for most of 
the cruciform region of acceptable relic density, and B s — > n + fi~ has no significant impact. 

These two examples illustrate that the extra degree of freedom in the NUHM2 opens up 
even broader possibilities for reconciling supersymmetry with the LHC constraints. 

5 Results for sub-GUT CMSSM models 

We now consider a different one-parameter extension, one in which the soft supersymmetry- 
breaking parameters mi/ 2 , m and A are assumed to be universal at some input scale M in 
below the conventional supersymmetric GUT scale, as might arise in models where the 
supersymmetry-breaking dynamics acts below the GUT scale. As was discussed in [21 1 and 
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Figure 8: The NUHM2 (p, m A ) planes for m 1/2 = 1000 GeV, m = 1000 GeV, A = 2.5 m , 
and tan /3 = 10 (left) and tan ft = 30 (right). The interpretations of the shading and contour 
colours are described in the text. 



summarized earlier, the sparticle spectrum in such a scenario is in general more compressed 
than in the conventional CMSSM with, e.g., smaller mass differences between squarks and 
sleptons and between different types of inos. Models with similarly compressed spectra and 
the prospects for their discovery have been studied in [42] . 

As presented in [21], this compression of the spectrum causes the relic cold dark matter 
density to decrease at generic fixed values of (m 1 / 2 , toq) as M in decreases, and hence the strips 
of parameter space with the appropriate dark matter density tend to move inwards from the 
boundary of the region allowed by the neutral LSP and electroweak symmetry breaking 
constraints. These features are visible in the upper panels of Fig. |9j which are (m^^o) 
planes for tan/3 = 10 and 40, both with M in = 10 11 GeV, A Q = 2.5 m , and /i > 0. In both 
cases there are wedges of the (mi/2,m ) plane at small and large m (and hence large 
A ) that are excluded because the lighter stop, t\, is the LSP in addition to the common 
wedge in the lower right of the planes where there is a f\ LSP. We see that in both panels 
there is a prominent funnel due to s-channel annihilations of LSPs through the heavy Higgs 
scalar and pseudoscalar, where the relic cold dark matter density is suitable. There are 
also other very thin strips running above and roughly parallel to the prominent funnels, due 
to rapid s-channel coannihilation processes between the LSP and heavier neutralinos and 
charginos. In both cases, we also see t\ coannihilation strips near the upper boundary of the 
allowed \ LSP region, in the neighbourhood where focus-point strips appear in the CMSSM. 

We now discuss how these sub-GUT models fare with the available LHC constraints. The 
MET constraint (not shown) excludes only a small corner of the displayed triangle between 
the fx and t\ LSP boundaries. As usual, B s — > / u + ^ does not impact the allowed parameter 
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Figure 9: Sub-GUT CMSSM scenarios with tan/3 = 10 (left) and tan/3 = 40 (right) in 
both cases with A = 2.5m , and for M in = 10 11 GeV (upper) and 10 9 GeV (lower). Plotted 
contours of Mh are 114, 119, 122.5, 124, 125, 126, and 121 GeV. The interpretations of the 
shading and contour colours are described in the text. 



space for tan/3 = 10, but does have significant impact for tan/3 = 40, requiring mx/ 2 ^ 1500 
to 2000 GeV along the prominent rapid-annihilation funnel, and mi/s ^ 1000 to 1500 GeV 
along the stop coannihilation strip. The most novel feature of these sub-GUT models is 
the relative ease with which they respect the Higgs mass constraint. Most of the funnel for 
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tan (3 = 10 is compatible with Mh = 125 to 126 GeV, once the uncertainties in the theoretical 
calculation within the CMSSM are taken into account. This is also true along the funnel 
when tan/3 = 40, for the points compatible with B s — > On the other hand, in both 

cases the value of Mh is too small along the stop coannihilation strip. 

The lower panels of Fig. [9] display (mi/2,mo) planes for the same values of tan/3, /i > 
and Aq = 2.5 mo, but now with Mj n = 10 9 GeV. In each case, rotating anticlockwise around 
the plane from the lower right, we first see a stau coannihilation strip followed by a rapid- 
annihilation funnel (which are both more prominent for tan/3 = 10). In the tan/3 = 10 case 
(lower left) we also see 'echo' funnel strips where is smaller for the same value of m , 
which are due to s-channel coannihilations between the LSP and more massive neutralinos 
and charginos. At lower Mj„, there is less running and stops remain heavier than the lightest 
neutralino at lower mo, thus, we see a recession of the stop LSP region at lower Mj n . Next 
to this region, we see a prominent stop coannihilation strip, but note that only for small m 
does it have an acceptable value of Mh- In the tan /3 = 40 case, the s-channel funnels merge 
with the stop coannihilation strip for (mi/2,m ) ~ (1000,2000) GeV (adjacent to a region 
excluded by b — > 57), again yielding acceptable Mh for lower ttlq. The LHC MET constraint 
is not relevant in either of the lower planes, and neither is B s — > fi + when tan (3 = 10. 
However, B s — > n + fi~~ does exclude almost all the dark matter regions for tan /3 = 40, with 
the exception of the stau coannihilation strip which reemerges at very large m^j. 

These examples illustrate how sub-GUT CMSSM models are better able to accom- 
modate the LHC constraints, principally because the greater importances of additional 
(co) annihilation processes permit larger values of and mo than in the conventional 
GUT-scale CMSSM. 



6 Results for sub-GUT mSUGRA models 

Encouraged by the relative success of sub-GUT CMSSM models in accommodating the latest 
LHC constraints together with other constraints, we now consider whether this compatibility 
can be maintained if we remove one of the sub-GUT CMSSM parameters by imposing the 
mSUGRA relation between the trilinear and bilinear soft supersymmetry-breaking parame- 
ters: Aq = i?o + mo, and also the mSUGRA relation between the gravitino and scalar masses: 
m 3/2 — m o before renormalization. 

We start with the (mi/ 2 ,m ) planes for the Polonyi case A = (3 — \/3) m shown in 
Fig. 10 (Results for a GUT-scale Polonyi model were shown in the left panel of Figure [3]) 



The case M in = 10 13 GeV is displayed in the upper left panel of Fig. 10 , where we see that 
there is a forbidden region at large mi/2 and mo ~ 2000 GeV where the fi is the LSP. At 
larger and/or smaller m.1/2 the lightest neutralino \ is the LSP, and at smaller m the 
gravitino is the LSP. There is a narrow funnel in the \ LSP region with suitable cold dark 
matter density that barely reaches (mi/2,mo) ~ (2000, 1500) GeV, and a 'boomerang' region 
at lower mo, extending to large mi/2, where the density of gravitino LSPs produced in NLSP 
decays is suitable for the cold dark matter. The tip of the \ LSP funnel and the upper 
arm of the gravitino boomerang are just barely compatible with Mh ~ 124 GeV, which is 
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compatible with the LHC measurement within uncertainties. We note that all (most of) 
the visible part of the plane is compatible with the LHC B s — > /j + //~ (MET) constraint, 
including these two M/j-compatible regions. In this plane, tan (3 is relatively small except 
when mi/2 is very large. This accounts for the sufficiently low value of B s — > fi + fi~ relative 
to the Standard Model (between 1.0 and 1.3 over much of the displayed plane). 

Turning now to the upper right panel of Fig. 10 for Mj n = 10 11 ' 5 GeV, and descending 
from the boundary of the electroweak symmetry breaking region in the upper left to lower 
right, across the \ LSP region, we encounter several new dark matter-compatible strips, 
some of which were present embryonically in previous figures. The uppermost in a series of 
wisps is due to X3 ~ X ± coannihilation enhanced by a near-on-shell t-channel pole. Moving 
across the plane towards larger 777.1/2, several wisps are crossed prior to reaching the most 
prominent diagonal feature, which is the more familiar x~ X s-channel funnel. These wispy 
features above the x ~ X funnel are due to coannihilations of neutralinos and charginos 
enhanced by s-channel H/A poles. Further down and to the right there is a coannihilation 
strip adjacent to the forbidden fi LSP triangle. Continuing below this triangle, we encounter 
a near-horizontal band where fi decays yield a suitable density of gravitino LSPs. In this 
case, almost all the (m^j, mo) plane is compatible with the LHC MET and B s — > 
constraints, and large parts of the many dark matter-compatible strips in the x LSP region 
are compatible with the measurement. 

To help understand the origin of the multiple avenues for coannihilation, we show in 
the upper left panel of Fig. fTTI the sparticle spectrum as a function of mi/ 2 for fixed m 
= 2000 GeV for Mj„ = 10 GeV. Here we see clearly the possibility of multiple resonant 
reactions when m^/2 is equal to m X3 at 777-1/2 ~ 800 GeV, m X2 at 777.1/2 ~ H00 GeV, and 
m x at 777-1/2 ~ 1800 GeV. In the upper right panel of Fig. 11, we see the behaviour of the 
neutralino composition as a function of the universal gaugino mass: the Higgsino composition 
is enhanced at large mi/2, contributing to an increased cross section and hence acceptable 
relic density past the point where the stau coannihilation strip would normally end. 

Looking back at the (777,1/2, m o) plane for M in = 10 ia5 GeV shown in the lower left panel of 
Fig. 10, and again descending across the x LSP region away from the electroweak symmetry 
breaking boundary towards the f\ LSP region, we again encounter several near-parallel dark 
matter-compatible strips. The first is a relatively broad focus-point region where the LSP 
is more Higgsino-like, then there is a coannihilation funnel, followed by the more familiar 
X — X rapid- annihilation funnel. Finally, we encounter a Higgsino - f\ coannihilation strip. 
Below the f\ LSP region, there is again a band where fi NLSP decays to gravitinos yield 
a suitable relic dark matter density. As in the previous case, most of the (777-1/2, m o) plane 
is again compatible with the LHC MET and B s — > constraints due to the relatively 

large mass scales and/or small value of tan/3. The upper parts (mo 2 TeV) of the dark 
matter-compatible strips in the x LSP region are compatible with the Mh measurement. 

Finally, we consider the (mi/ 2 ,m ) plane for M in = 10 10 GeV shown in the lower right 



panel of Fig. [10] At large mo, we see that the broad focus-point strip is joined to one side 
of the standard x ~ X funnel. The other side of this funnel is connected to a Higgsino 
coannihilation strip hovering above the f\ LSP boundary. At lower mo and mi/2, there is 
another loop connecting the lower part of the focus-point strip with the X~X funnel, and at 
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Figure 10: The (mi/2,mo) planes for fi > in the sub-GUT Polonyi mSUGRA model with 
A = (3 - \/3)m and M in = 10 13 GeF fupper left), 10 1L5 GeF fupper r^/i^ ; 10 ia5 GeV 
(lower left) and 10 10 GeV (lower right), found using the latest version of the SSARD code fSty . 
The interpretations of the shadings and contour colours are described in the text. 



low mo but larger m\/2 there is again a band where gravitinos could provide the dark matter. 
As before, most of the (roi/2, mo) plane is compatible with the LHC MET and B s — > /j + fi~ 
constraints, and there are large stretches of the dark matter strips at large and mo that 
are compatible with the LHC Higgs measurement. 
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Figure 11: The sparticle spectrum (left panels) and (right panels) the composition of the 
LSP x functions ofm 1 / 2 in a sub-GUT Polonyi mSUGRA scenario with m = 2000 GeV, 
for Mi n = 10 11 ' 5 GeV (upper panels) and = 10 10 GeV (lower panels). 



For this case also, we display the sparticle masses and neutralino composition as functions 



of mi/2 for mo = 2000 GeV in the lower panels of Figure 11 Here, we clearly see the 
congruence of masses at 7711/2 ~ 1000 GeV all contributing to a resulting low relic density. 
In this case, we also see that, for this value of mo, the neutralino is predominantly a Higgsino, 
but in the funnel centre, there is an enhanced bino component leading to a very mixed LSP 
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which also leads in general to an enhanced cross section. 

The above discussion was for sub-GUT Polonyi mSUGRA models, and we now discuss 
briefly what happens in sub-GUT mSUGRA models with values of A ^ (3 — \/3) m . First, 
we recall that in mSUGRA models the value of tan is fixed by the electroweak conditions 
as a function of the other, independent model parameters, as shown by the grey contours 
in the different panels of Fig. 10 As shown there, in the Polonyi case typical values are 
tan/3 ~ 20 to 30, which is why the B s — > fi + n~ has relatively little impact - though this 
does increase slightly at smaller M in where larger values of tan /3 are found. On the other 
hand, we have already seen in the right panel of Fig. [3] that the derived values of tan /3 are 
as large as ~ 40 for M in = Mqut- 

In sub-GUT mSUGRA models with Aq = 2 mo, as shown in Fig. 12 we see that tan/3 
increases to ~ 45 to 50 for Mj n = 10 14 and 10 13 GeV. This has two consequences worth noting: 
(i) the impact of the B s — > / u + pT constraint is increased to such an extent, compared to the 
Polonyi case, that no A = 2m sub-GUT mSUGRA models are valid for either value of 
M in , and (ii) for M in < 10 13 GeV, the increasing value of tan/3 induces non-perturbative 
instabilities in the renormalization-group equations over much of the (m^jWo) plane. For 
these reasons, sub-GUT mSUGRA models with Aq much larger than the Polonyi value are 
disfavoured. 
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Figure 12: The sub-GUT mSUGRA (mi/2,mo) planes for A = 2m and M in = 10 GeV 
(left), Mi n = 10 13 GeV (right). The interpretations of the shading and contour colours are 
described in the text. 



In the case of sub-GUT mSUGRA models with Aq less than the Polonyi value, we recall 
(see also the left panel of Fig. [3j, that Mh is much smaller than the LHC measurement if 
A = and M in = Mqut- This problem persists for M in < Mqut, so we do not discuss 
further this possibility. 
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The key messages from this study are that (i) many (co) annihilation processes that are 
relatively unimportant in the CMSSM can become important in sub-GUT models, as the 
sparticle spectrum becomes more degenerate and Boltzmann suppression factors therefore 
become less significant. As a result, (ii) the cold dark matter density may fall into the 
favoured range almost anywhere in the interior of the (mi/ 2 ,m ) plane away from the f\ 
LSP and electroweak symmetry breaking boundaries, depending on the precise value of M in 
and A . Next, (iii) many of the regions allowed by the dark matter constraint are compatible 
with the current LHC constraints, including the measurement. Finally and intriguingly, 
(iv) in the specific cases of sub-GUT mSUGRA models, values of A very different from the 
Polonyi value A = (3 — V3) m are disfavoured. 

7 Summary 

It is well known that only very small regions of the CMSSM parameter space are consistent 
with all the theoretical, phenomenological, experimental, astrophysical and cosmological 
constraints, notably including the preferred range of the cold dark matter density. We have 
explored in this paper the extent to which some generalizations of the CMSSM may wriggle 
out of the straitjacket imposed, in particular, by the LHC constraint on missing-energy 
(MET) events, by the LHCb measurement of the branching ratio for B s — > fi + fi~ and the 
LHC measurement of the mass of the Higgs boson. The principal classes of models studied 
are those with non-universal supersymmetry-breaking contributions to the Higgs masses 
(NUHM1,2) and models in which universality is imposed on the supersymmetry-breaking 
parameters at some input scale M in < M GUT (sub-GUT CMSSM and mSUGRA). 

It is not difficult to find regions of the NUHM1,2 parameter spaces where the cosmological 
cold dark matter density falls within the preferred range, even if the sparticle masses are 
relatively large, as required by the LHC MET and Mh constraints. For example, this may 
happen in a transition region where the x LSP has a relatively large Higgsino component. In 
the CMSSM, this could happen only in the focus-point region, a possibility now disfavoured 
by the XENON100 upper limit on cold dark matter scattering, but in the NUHM1 the \i 
parameter may be regarded as free and the transition region may appear in other parts of 
the (mi/2,iTio) plane where the LHC MET and Mh constraints are respected and BK(B S — > 
is close to the Standard Model value. 

Sub-GUT models offer other possibilities for reconciling the LHC and other constraints. 
The compression of the spectrum as M in decreases implies that more coannihilation pro- 
cesses may become important, suppressing the relic density below the range expected in the 
CMSSM, and enabling it to be compatible with the cosmological constraint also at relatively 
large values of mi/ 2 and m . Also, in sub-GUT models we find a proliferation of strips where 
the dark matter density falls within the preferred range. 

The NUHM1,2 and sub-GUT models are examples to illustrate that there is no general 
conflict between the LHC constraints and supersymmetry, even within the MSSM frame- 
work. It suffices to relax the constraint of universality of the soft supersymmetry-breaking 
parameters at the GUT scale, as we have demonstrated in a few types of scenarios. How- 
ever, it is worth noting that the surviving models in these particular scenarios tend to have 
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quite large sparticle masses, and do not provide a supersymmetric resolution of the — 2 
discrepancy. 
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